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Abstract 

We compute the strong coupling constant Gb** Bit (Gd** d-k), where B** (D**) 
is the + P— wave bq (cq) state, by QCD sum rules and by light-cone sum 
rules. The two methods give compatible results in the limit mq — > oo, with 
a rather large value of the coupling constant. We apply the results to the 
calculation of the hadronic widths of the positive parity B and D states and 
to the chiral loop contribution to the ratio /d s / Id- 
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I. INTRODUCTION 



In the light quark (q = u, d, s) zero mass limit (m q — > 0) and in the heavy quark (Q = c, b) 
infinite mass limit (tuq — > oo) Quantum Chromodynamics exhibits symmetries that are not 
present in the finite mass theory: chiral SU(3)l x SU(3)r symmetry, heavy quark spin and 
flavour symmetries M, as well as the velocity superselection rule ||, valid insofar only strong 
interactions are considered. 

All these symmetries can be used to build up an effective chiral lagrangian [H] for light 
pseudoscalar mesons and heavy (Qq) negative parity (J p = 0~, 1~) mesons. This effective 
lagrangian contains a symmetric term plus corrections to the heavy-light symmetries such 
as, for example, terms proportional to powers of 1/rriQ or terms proportional to powers of 
the light quark masses m q . In the spirit of the chiral effective theory, the resulting lagrangian 
is also an expansion in the light meson fields derivatives. 

Besides the octet of the pseudo-Goldstone bosons 7i,K,r] and the (D, D*), (B, B*) 
states, one can include the light vector mesons belonging to the low lying SU(3) nonet: 
p, K*, us, uj , using the so-called hidden symmetry approach [§]||. As well known, in this 
approach || the lagrangian exhibits an extra local gauge symmetry SU (3)#, and the 1~ light 
meson octet represents its gauge bosons. They acquire a mass because SU{3)h is sponta- 
neously broken. Quite recently it has been observed [0] that, if SU(3)h is unbroken, a new 
symmetry (vector symmetry) arises. Its implications for the heavy-light chiral lagrangian 
have been examined in Ref. f|. Additional notions arise, in general, when also axial- vector 
bosons are present [f|. 

Another extension of the heavy-light lagrangian is obtained by including effective fields 
describing positive parity (Qq) mesons. According to the value of the angular momentum of 
the light degrees of freedom, (sf = | + , | + ) the Heavy Quark Effective Theory pre- 
dicts the existence of two multiplets, the first one comprising + and 1 + mesons, the second 
one containing 1 + , 2 + states. The role of the \ + doublet (0 + , 1 + ) in some applications of 
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chiral perturbation theory has been considered in [|T^][|. Another application is in the realm 
of the semileptonic D and B decays J3|. 

The aim of the present paper is to give an estimate of the strong coupling constant 
describing the interaction of the pseudoscalar light mesons with the positive parity (sf = 
| + ) and negative parity (sf = | ) heavy mesons and to make quantitative estimates of the 
effect of the positive parity heavy mesons in some calculations in chiral perturbation theory. 

After a review of the heavy-light chiral lagrangian in Section II, we consider two sum 
rules for these coupling constants: the first one, based on the method of the single Borel 
transform (in the soft pion limit), is discussed in Section III, while in Section IV we derive 
these couplings by the method of the light-cone sum rules (for a review see [|14|j ). 

In Section V we compute the decay widths of the excited positive parity heavy mesons, 

1 ~f~ ^ H~ 

and we estimate the mixing angle between ^ and | axial vector states. 

In Section VI we comment on the role of the positive parity heavy mesons in the chiral 
loop contributions to the ratio /d s //d, which, as observed in [|12[], may be considerable. We 
find the numerical result: 

& = 1-09 (1.1) 
Id 

where part of the SU(3) violation effect /d s //d 7^ 1 m ay be attributed to positive parity 
states. 

Finally, in Section VII we draw our conclusions. 

II. THE HEAVY-LIGHT CHIRAL LAGRANGIAN 

In the effective heavy-light chiral lagrangian, the ground state (sf = | ) heavy mesons 
are described by the 4x4 Dirac matrix 



x There is no coupling of the (sf = § + ) states to the (0,1) heavy meson doublet and to the 
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pseudoscalar Goldstone bosons at the lowest order in the chiral expansion 
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[PLY - Pal,] 



(2.1) 



where v is the heavy meson velocity, P* M and P a are annihilation operators of the 1~ and 0~ 
Qq a mesons (a = 1, 2, 3 for u, d and s): for charm, they are D* and D respectively. Similarly, 
the positive parity 1 + and + (sf = | + ) are described by 

1+f 



S a 



(2.2) 



It should be observed that all the operators appearing in Eq. (|2.1| ) and Eq. (|2.2j ) have 
dimension | since they contain a factor ^mp in their definition. 

As for the octet of the pseudo Goldstone bosons, one uses the exponential form: 



where 



M 



2 



7T 



V 



K~ 



%M 

i = exp — 

JIT 



7] 7T + K + 



(2.3) 



(2.4) 



and U = 132 MeV. 

The lagrangian describing the fields H, S and £ and their interactions, under the hypoth- 
esis of chiral and spin-flavour symmetry and at the lowest order in light mesons derivatives 
is: 

C= f -f< mdptf >+i< H b v»D^ ba H a > 

o 

+ < S b (i v^D^ ba - 5 ba A)S a > + ig< H b ^ 5 A^ a H a > 

+ ig' < SiJ^A^Sa > + [lh< SbJ^A^Ha > + h.c] (2.5) 



where < . . . > means the trace, and 



ba 



(2.6) 
(2.7) 



E = £ 2 and A is the mass splitting of the S a states from the ground state H a . Numerically 



we use A = 500 ± 100 MeV, an estimate based on quark model |T5j and QCD sum rules 
|16| computations of the masses of the excited sf = | + mesons. 

In Ref. |T7| a QCD sum rule has been considered to compute the strong coupling constant 
g appearing in Eq.([2.5|), with the result (at the order a s = 0): 



# = 0.39 ±0.16. (2.8) 
This result, valid in the soft pion limit, has been confirmed by a subsequent, independent 



analysis, based on the method of light-cone sum rules |L8 



g = 0.32 ±0.02. (2.9) 

The results of similar sum rules for the strong coupling constant h in Eq.( [2.5|) will be 
presented in the subsequent Sections. 

III. QCD SUM RULE FOR h 

Let use define the strong amplitude 

G B ** Bn = < n+(q) B°{q 2 )\B**+( qi ) > (3.1) 

where B** is the + state in the sf = \ + doublet. The amplitude Gb**Bw is related to 
the strong coupling constant h appearing in the heavy-light chiral lagrangian (|2.5|) by the 
formula: 

m? B ,* — m 2 B h 

G B **Bw = - Jm B m B ** -j- ■ (3.2) 

m B ** f w 

In the limit — > oo one has: 

m B = m b + oo ± O (—) 
KrribJ 

m B ** -m B = A + O (— ) (3.3) 

\m b J 
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and 

2h 

Gb**Btt ^- T m b A (3.4) 

In order to derive a sum rule for Gb** b-k and h, we consider the correlator: 

A, = ijdx< n+(q)\T(j 5 (x)V,(0)\0 > e~^ x = A(ql q 2 2 , q 2 )q, + B (qf , q 2 2 , q 2 ) P, (3.5) 

where % = uij^b, = b^^d, q = q x — q 2 and P = q\ + g 2 . The scalar functions A and 
B satisfy dispersion relations (D.R.) that are computed, according to the method of QCD 
sum rules, in two ways: either by saturating the dispersion relation by physical hadronic 
states, or by means of the operator product expansion (O.P.E.). In Ref |T7| we considered 



the dispersion relation for the scalar function A and we used it to compute the coupling 
constant Gb*b-k defined by the matrix element: 

< 7r+(g) B°(q 2 )\B*+( qi ,e) > = G b *b. H ■ <f . (3.6) 

We obtained the result: 

f B fs* G B *B, = 0.56 ± 0.12 GeV 2 (3.7) 
while for charm we obtained: 

fo fn* G D * Dn = 0.34 ± 0.08 GeV 2 . (3.8) 
The leptonic decay constants appearing in the previous equations are defined as follows: 

< 0\h^ 5 d\B(p) > = if B jf 

< 0\wfb\B*fa e) > = i m B *f B * ^ . (3.9) 

It should be noticed that the phases in the previous equation are consistent with the defini- 
tion of the weak current in the effective theory, see below Eq. ( |6.3j ). 
In the following we shall make use also of the matrix element: 

< 0\uYb\B** (p) >= if B ~ p" . (3.10) 



It may be observed that the hadronic side of the sum rule for A also includes the B** pole in 
the variable q\ (the low lying pole in the variable q\ is provided, of course, by the B meson). 



In [|I7]] one gets rid of it by exploiting the tiny mass difference between B and B*, which 
allows to include the J p = + pole in the so called parasitic terms that may be appropriately 
parametrized. Here we are interested precisely in this pole and, in order to obtain it, we 
consider the dispersion relation for the scalar function B in Eq. fl3.5|) that should be written 
in general as follows fl9| : 

p(s, s',q 2 ) 



B{q\,ql,q 2 ) = ^ J dsds'- 



s-ql){s'-ql) + 

„2\ „ f n „2\ 



+ p i\q\) \ ds — — + P2{q 2 ) / ds — + P 3 {q 1 ,q 2 , q 

J S — Qn J S — Of 



,2 „2 „2^ 

- ... , " r --sl< ' 

q$ J s- q{ 

(3.11) 



As proven in p0[ Pi{q\) = -fa (^i) = ^ since only this value of the polynomials Pi(qf) 
and ^2(02) i s compatible with the vanishing of the form factors for large values of their 
arguments, as predicted by quark counting rules. The subtraction polynomial P^q^.q 2 ) 
does not contribute to the sum rule since it vanishes after the Borel transform; therefore we 
shall neglect it in the sequel. 

In order to compute the D.R. for B, in this Section we make the approximation of the 
soft pion limit (S.P.L.): q — > 0. This approximation presents the advantage of a consider- 
able simplification of the calculations; moreover in this scheme the limit itlq — > 00 can be 
performed in a well defined way. On the other hand it might be argued that in the decay 
B** — > Bit the pion momentum is not constrained to be small, since m#** -mg ~ A ~ 500 
MeV. We shall comment on the uncertainties introduced by a small pion momentum ap- 
proximation in the next Section, where we shall present a light-cone sum rules calculation 
of Gb**btt which is not based on S.P.L. 

The problem arising in the soft pion limit is related to the fact that q = implies 
q\ = q\. As a consequence, one cannot perform a double Borel transform in the variables 



qf , gl, and the single Borelization procedure has to be used. As well known |TU] in this way 



one introduces unwanted not-exponentially suppressed contributions (the so called parasitic 



terms) that have to be estimated^. We shall show in the sequel how this problem can be 
solved. For the time being we consider the result of computing B(qf, qf, 0) by OPE in the 
soft pion limit. The result of a straightforward analysis (much similar to that considered in 



Refs. [T^J2T[) is as follows: 

B{ql qf, 0) = B(°) + B« + + B^ + + B® 



(3.12) 



where 



B {A 

B {5 



< uu > 



fn 9l ~ ™\ 



TOq < UU > 





m\ <uu> 
±fMl ~ m 2 b ) 2 
. 



2\2 



1 - 



2ml 



q{-m% 



(3.13) 



In Eqs. (|3.13| ) < uu > is the quark condensate (< uu >= — (240MeV r ) 3 ), ml is defined by 
the equation 

< ug s o ■ Gu >= ml <uu> (3-14) 



and its numerical values is: ml = 0.8 GeV 2 . The origin of the different terms in Eq.( p.l2|) 
is as follows. B^ is the leading term in the short distance expansion; B^, B® and B® 
arise from the expansion of j$(x) at the first, second and third order in powers of x\ B^ 
and B^ arise from the expansion of the heavy quark propagator at the second order and 
from the zeroth and first term in the expansion of j$(x) respectively. We have considered 
all the operators with dimension D < 5 in the O.P.E. of the currents appearing in Eq. (|3.5| ). 

Let us now compute the hadronic side of the sum rule, that we call B^. We divide the 
integration region into two regions Di,D 2 , as depicted in Fig. 1. D\ is bounded by the lines 
s — ml, s' — ml and s + s' = C. We assume that C satisfies the bounds: 



2 A similar situation is met computing gB*Bn by this method [17]. 
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min(s'o + m B **, s + m 2 B ) >C> m 2 B „ + m\ (3.15) 

where So and s'q are thresholds for continuum production in the variables s and s' respec- 
tively. The bound in ( |3.15| ) is chosen in such a way that, inside Di, only the poles B* and 



B** (in the variable s) and B (in the variable s') are included. Their contribution is as 
follows: 

g ( 2 2 q\ _ G B .BjB/ffm fi (mg, — m\) G B ** Bn f B f B **m B /g 

P ° le u v 4rn b m B *(qf-rn B *)(q 2 -m B ) 2m b (q 2 - m 2 B *,){q 2 - m%) ' 

We observe that, due to the presence of the factor m 2 B * —m 2 B the contribution of the 1~ pole 
is strongly depressed as compared to the term describing the + resonance; in the rrib — > oo 
limit, the ratio of the two terms is 0{-^). We have checked that this suppression holds not 

b 

only for the beauty, but also for the charm, where the vector meson pole contribution is less 
than 3% of the + pole. For this reason we shall omit the 1~ pole in the sequel and we shall 
take only the scalar contribution in Eq. (|3.16 ). 



Let us now consider the region D2, where a continuum of resonances contributes. Fol- 



lowing Ref. |2z] we assume the following model for the hadronic dispersive function in this 
region: 

p(s,s') = f(s,s')8(s- s')9(s + s' - C) [region D 2 ) . (3.17) 

In order to justify this choice let us introduce the variables z = s — s' and y = In these 
new variables we write the hadronic side of the sum rule as follows: 

BhaM, <l!, 0) = B pde (q 2 , qf, 0) + B cont [q\) (3.18) 

where the contribution of the continuum of resonances (region D2) is as follows: 

1 f°° / > +2(j/— m?) I Z Z 

B cont {ql) = — / , d y / , 9N dz—— z 2\P(V+ o»V- o) • ( 3 - 19 ) 

■n 1 Jc/2 J-2( y -m 2 b ) [y + § - qt){y - f - qf) 2 2 

Let us now rewrite the leading term in the O.P.E. expansion as an integral over the variables 

z, y- 
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U Qi ~ m b 



< uu > r°° 



dy- 



+2(y-m 2 b ) 



dz8(z) . 



[3.20) 



To relate (|3.19|) and ( [3.20D , one imposes local duality in y [22]; this means that inside the 



region D 2 , i.e. for y > C/2, the arguments of the ^-integrals in Eqs. ( 3.19 ) and ( 3.20 ) 
should be equal. This implies that ( p. 17 ) is valid with 



f(s,s') 



, < uu > 



71 



fn 



(3.21) 



which fixes the model of the hadronic continuum in the region D 2 . A justification for local 



duality in the variable y can be obtained by explicit calculations in particular models [24 . 

Using these results we obtain, after the Borel transform, the following sum rule, which 
as we have already observed, is valid only in the soft pion limit: 

G B** Bit f B f B** m B 



m B * t — m B 



-m? B /M 2 



-m 2 B ,JM 2 



2ml, < uu > 



e -m 2 /M 2 



1 1 

AM 4 



C/{2M 2 ) 



(3.22) 



The sum rule is valid provided one imposes two conditions on the Borel parameter M 2 . First 
one requires that the contribution of the continuum does not exceede the pole contributions 
(for larger values this would produce uncontrollable uncertainties in Gb**b-k due to our poor 
knowledge of the continuum); this fixes the upper bound for M 2 . On the other hand, for 
the O.P.E to be meaningful, we have to impose that the higher power terms in 1/M 2 in 
Eq.( p.22|) have decreasing values, which fixes the lower bound for M 2 . For the beauty sector 
we use the values m#« — m# = A pa 500 ± 100 MeV, m& = 4.6 GeV, the parameter C in 
the range (61 — 64) GeV 2 (with thresholds s = s' ~ 36 GeV 2 ). The previous criteria are 
satisfied for M 2 in the range (15 — 25) GeV 2 . We get the result: 



G b** Btt f b f b** = 0.43 ± 0.06 GeV 3 . 



(3.23) 
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For charm, we use mo** — mp & A = 500 ±100 MeV, m c = 1.35 GeV, and the parameter 
C in the range(10 — 14) GeV 2 (with thresholds Sq = s' — 6 GeV 2 ). The criteria for M 2 are 
satisfied in the range M 2 = (2 — 8) GeV 2 , and one gets the result: 



Gb-dJd/b.. = 0.38 ± 0.11 GeV 3 . 



(3.24) 



To get G B **Bn and G D ** Dn we use f B = 0.18 ± 0.03 GeV and / B « = 0.18 ± 0.03 Gel/ 
16fl obtaining: 



Gb** Bit = 13.3 ±4.8 GeV 



(3.25) 



while, using f D = 0.195 ± 0.020 GeV and fo, = 0.17 ± 0.02 GeV, JT6[ we obtain: 



Gd** d-k = 11.5 ±4.0 GeV . 



(3.26) 



We observe substantial violations of the scaling law Gd**d-k/Gb**b-k ~ On the other 
hand the ratio i? = g B ** g7r j^j B ** is less sensitive to scaling violations (numerically we find 
R ~ 1.13, to be compared to the scaling prediction R — 1). 

Let us now take the limit m& — > oo. This limit is defined by Eqs. ( |3.3| ) and by a rescaling 
of the Borel parameter M 2 : 



M 



2m h E 



(3.27) 



In this way we obtain the asymptotic rule: 



hFF~ 



-A/E 



4 <UU> 



3 W/E (1 



-6/E 



16E 2 ' 



(3.28) 



the quantities 5, F and F + , defined by: 



C 

m 



B 



F = fBy/™B 



F + = f B **V^ 



(3.29) 



remain finite in the infinite heavy quark mass limit, modulo logarithmic corrections. The 
constraints on the Borel parameter imply that E must be in the range (0.5 — 1.4) GeV . 
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Using u = 0.50 ± 0.07 GeV, A = 500 MeV, and 5 = 400 ± 50 MeV the numerical outcome 
of the sum rule is: 



hFF A 



0.072 ± 0.008 GeV z 



(3.30) 



(the results are weakly dependent on 5). By the values F = 0.30 ± 0.05 GeV 3//2 |?3 and 
F + = 0.46 ± 0.06 GeV 3 / 2 @ we obtain: 



h = -0.52 ±0.17 . 



(3.31) 



IV. A LIGHT-CONE SUM RULES CALCULATION OF h 

An independent calculation of the strong coupling h can be carried out using a method 
that allows us to consider arbitrary momenta of the pion. We shall consider the correlation 
function 

A(qlqlq 2 ) = i j dx < n + (q) \T(j 5 (x)j(0) |0 > e"^ (4.1) 

where j'5 = iuj^d (as before) and j = bd. 

The method consists in expanding the T-product of the quark currents, appearing in 
Eq.( [4.1|) , near the light-cone, in terms of non-local operators whose matrix elements can be 
written as wave functions of increasing twist. 

This approach finds its origin in the analysis of hard exclusive processes in QCD ||13|| : 
within this framework, strong couplings (g^pir, QttNn), form factors (such as tcA^*, the pion 
form factor at intermediate momentum transferred) and nucleon magnetic moments have 



been calculated . Light-cone sum rules have also been used to calculate the form factors 



governing B and D meson semileptonic decays ||25|| , and the radiative B —>■ K*^ transition 



In our calculation of h we follow the notations adopted by the recent paper by Belyaev 
et al. [TH| devoted to the calculation of the coupling constants gB*Bn and gD*Dw Within the 



light-cone sum rules approach the correlator Eq. (|4.1|) can be written as follows: 
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Jo ^ ml - (q 2 + uq) 2 [ml - (q 2 + uq) 2 } 2 ^ m\ - {q 2 + uq) 2 > > 

+ M?2 -q) du Va iT - — ■ -— - . 4.2 

Jo J [mf - (q 2 + q(a x + ua 3 )) 2 \ 2 

The pion wave functions (pp(u), (p n (u) and g\{u) appear in the matrix elements of non- 
local quark operators p~8f: 



< 7r(g)|J(a;)7 M 75«(0)|0 > = -if^ f due^&M + x 2 9l (u) + 0{x A )) 

J 



2 ^ 

+ UK ~ ^zr) [_ due^ x g 2 {u) , (4.3) 



qx Jo 



< 7i{q)\d{x)n 5 u{0)\0 > = Jn g / dMe m "> P (M), (4.4) 

where is related to the light-quark condensate < uu > by the current algebra 

relation ;7rm,n = —2 <u t u> . It should be noticed that the path-ordered gauge factor 

m u +m d f„ r a a 

Pexp (ig s Jq 1 dux^ A^(ux)) has been omitted in the above matrix elements due to the choice 
of the light-cone gauge x^A^x) = 0. 

The functions x{ a i^ a 2^z) and x(ai,a 2 ,a 3 ) are combinations of twist-four wave func- 
tions <f(ai, a 2 , a 3 ) and <p(ai, a 2 , a 3 ) 



x(«i, a 2 , a 3 ) = 2tp_ L (a 1 , a 2 , a 3 ) - f\\(ai, a 2 , a 3 ) (4.5) 

and 

x(«i, a 2 , a 3 ) = 2<p±(ai, a 2 , a 3 ) - tp\\(ati, a 2 , a 3 ) , (4.6) 
where the functions (p and (p parametrize the matrix elements of quark-gluon operators: 
< 7r(g)|d(x)7 M 755f s G a/3 (w£)w(0)|0 >= 

+U J ^(q a x p - qpx a ) I T>a m (a i )e i ^ + ^ (4.7) 
q ■ x J 

and 
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< n(q)\d(x)j^g s G al3 {ux)u(0)\0 > 



Qfi\9an 

q ■ x 



XaQ.fi 

q ■ x 



q ■ x 



D« f ^ i (a 1 )e i?l(£,1+ " 3) 



(q a x 3 - qpXa) / T>am\(ai)e 



iqx(ai+uct3) 



(4.8) 



In Eq. ( ^8|) the field G a p is the dual of G a p: G a p = ^e a 0s P G Sp ; the integration on the 
variables «j is performed considering that T>oti = da\da2da^5{l — a± — a 2 — 0:3). 

Using the normalization of the function ipp: Jq (pp(u)du = 1 we recover from Eq. Q4.2|) , 
in the soft pion limit q — > 0, the leading order expansion of B(qf, q 2 , 0) in Eq. (|3.13| ), apart 
from an overall factor —2m b due to the the different choice of the current interpolating the 
B** in the two cases (vector and scalar, respectively); it should be noticed that the terms 
proportional to the mixed quark-gluon condensate appearing in Eq. ( |3.13| ) are missing here 
since they are related to higher-twist pion wave functions. 

The hadronic representation of A(qf, yf , q 2 ) Eq. ( f4.1| ) can be expressed in terms of the 
contribution of the lowest-lying resonances B and B**\ 

1 



A P oie(qi,q2i ( l ) — Cb^btt/s/b* 



m 2 B m 2 B , 



(4.9) 



(m b + m d )(m b - m u ) (m|** - yf)(m| - yf ) 
in the region m 2 < s < sq, ml < s' < s' , and of the contribution of higher states and of 
the hadronic continuum. For q\ 7^ q 2 we can perform a double independent borelization 
in the variables —q\ and —q 2 . In this way, the parasitic contributions coming from the 
resonance-continuum terms in fl4.1|) are exponentially suppressed. From Eq. (|4.1| ) we get: 



B A 



had 



J ds ds' Phad (s, s>, q 2 ) ex P (-^ - JL) 



(4.10) 



where M\ is the Borel parameter associated to the variable — q\ and M 2 to —q\. In this 



expression the contribution of the pole reads 
1 



B A 



pole 



MfMf 



Gb^B-kIbIb* 



m 2 B m 2 B , 



(m b + m d )(m b - m u ] 



exp 



m 



B* 



(4.11) 



On the other hand, borelization of A® CD provides us with the following expression: 



BA QCD 



Ml Ml 



cxp 



M 2 



{- 2 



m b < uu > 



MVp(mo) 
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M 4 

- U— ti(u ) + 2U(M 2 + m 2 b )g[(u ) 



M 2 

~Y 

1 da 3 



u ° das 

[x( u o - Qt 3 , 1 - m , a 3 ) - x(uo - a 3 , 1 - u , a 3 ) 

o a 3 

(x(«o, 1 - Mo - «3, "3) + x(u , 1 - m - "3, «3)) 
1_a i (ia 3 



/o a 3 
2 / rfai , 

JO Juo—ai Q; 3 



x(«i, 1 - ax - «3, "3) } 



(4.12) 



where M 2 



and Mo = M 'i+ M -i , V 9 ' = ^ an d Pi = ^7- It is worth observing that, 
choosing the symmetric point Mo = \ (which corresponds to a quark and an antiquark of the 
same momentum inside the pion) Eq. ( |4.12|) is considerably simple, since at this point both 

(ft and g[ vanish. Moreover, at Mo = \ the subtraction of the continuum contribution can be 

2 2 

e - "^, at least in the twist 3 contribution WE 



we use 



done by substituting e - ^ 7 — > e - "^ 
this substitution everywhere in Eq. Q4.12 ) since higher twist contributions are numerically 
small). Therefore, we derive the sum rule: 



I exp 



2 2 

m B m B * 



exp 



[ 



m B *, + m 2 B 



(m b + m d )(m b - m u ) L 2M 2 
M 2 /^ r Am b < mm > 



- m£- 




so 1 


L M 2 - 


— exp 


- M 2 ^ 



(p P (u ) 



+ 



"0 C?Q/ 3 

[x(m - ar 3 , 1 - Uo, a 3 ) - x(«o - "3, 1 - u , "3)] 

a 3 

— (x(«o, 1 - uo - "3, "3)) + x(«o, 1 - M - a 3 , "3))) 
/o a 3 

/""° , ^«3 / si 

+ 2 y dax J — 2-x(ai, 1 - ai - a 3 , a 3 J j . 



(4.13) 



'«o— ai tt 3 

Let us now discuss the numerical analysis of Eq. Q4.13D - In the 6-channel we choose 
the same values for m b , rriB** and < uu > used in Section III. The main nonperturbative 
quantities are the twist 3 function (pp and the combinations \ and \ of twist 4 wave-functions. 
We choose the model in Ref . [27] , where such quantities have been fixed in the framework of 
a systematic expansion in the conformal spin. It turns out that, using y>p(l/2)| M6 = 1.07 [27 



the higher twist contribution in Eq.( |4.13| ) is 2% of the twist 3 contribution, and therefore 
the numerical value of <pp(l/2) represents a crucial quantity in our analysis. We allow the 
effective threshold s to vary in the range 36 — 40 GeV 2 . Moreover, we fix the highest 
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value of the Borel parameter M 2 in the duality window by imposing that the contribution 
of the continuum is 30% of the resonance. In this way we find M^ ax = 12 — 14 GeV 2 . The 
minimum value of M 2 is usually fixed by imposing that terms proportional to higher powers 
of 1/M 2 are small enough. Since in Eq.( |4.13p such terms are absent, we only look for a 
stability region in M 2 , and choose M^ in = 6 GeV 2 , which is the same value adopted in the 
analysis of gB*B-K- With these input parameters we find: 

G b** Bit f b f b** = 0.69 ± 0.14 GeV 3 (4.14) 

where the uncertainty is due to the variation of so and to the dependence of the numerical 
results on M 2 . Therefore, using the same values of the leptonic constants /# and /#»» 
adopted in the previous Section, we get: Gb**b-k = 21 ± 7 GeV and, from (|3.2j ), h{m b ) = 
-0.52 ±0.18. 

In the case of the charm sector, we use <^p(l/2)| Mc = 1.14, sq = 9 — 11 GeV 2 and M 2 in 
the range M 2 = 2 — 5 GeV 2 . We obtain: 

G D ** D J D f D ** = 0.21 ± 0.02 GeV 3 (4.15) 

and the results: Gd**d-k = 6.3 ± 1.2 GeV (using f D and jo** as in Section III) and h{m c ) = 
-0.44 ±0.09. 

A two-parameter fit of the above results in the form 

him) = h (1 + — ) (4.16) 
m 

gives the result: 

h = -0.56 ± 0.28 (4.17) 

and for the parameter a, a — 0.4 ± 0.8 GeV. 

Let us now compare these results with those of the previous Section. The values of h 
found by the two methods agree with each other. As for the finite mass results, the two 
methods sensibly differ (almost a factor of 2) in the case of the charm, while the deviation is 
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less important for the case of beauty (around 40%). These differences should be attributed 
to corrections to the soft pion limit that have not been incorporated in the results of Section 
III. 



V. EXCITED HEAVY MESONS WIDTHS 

In this Section we apply our results for the strong coupling constants to the calculation 
of the hadronic widths of the excited mesons. First of all, we can compute the strong widths 
of the decays Pq — > Ptt and Pi — > P*n. Differently from the decays of the positive parity 
states having = | + , i.e. the states P2 and P{, where the final pion is in P-wave, these are 



.S-wave decays; moreover, as shown in ||11|| , the single pion channels are expected to saturate 
the total widths. 

In the rriQ — > 00 limit one would obtain 

T(P - P+tT) = T(P 1 - P*+n-) = i- (jj A 3 (5.1) 

but this formula is of limited significance, especially for the case of charm, due to the large 
1 / itlq corrections coming from the kinematical factors. 
Keeping rriQ finite, the formulae become: 



1 



\'2\ 



[m 2 PQ - (ra P + m„) 2 )(mp o - (m P - mj 
Using G D * Dn = 6.3 ± 1.2 GeV, G B * BlT = 21 ± 7 GeV, and A D = A B = 500 MeV, one finds 



r ( p„ p+o = ^.. p ,^_^ - . (5.2) 



T(D -»■ Dn) ~ 180 MeV (5.3) 
T(B -> Ptt) ~ 360 Mel/ . (5.4) 

For the decay 1 + — >■ l~n (Pi — »■ P*7r) we find 

(m^ - (m P * + m 7V ) 2 )(mp i - (m P * - m^) 2 ) 



r(Pi -> p* + tt 



r 12 



87T 2771^ 



X 



1 / (m 2 Pi +mL 2 \ , . 

3 \ 4m 2 Pi m 2 P , J v 7 
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In the limit mg — > oo Eqs. fl5.5|) and (|5.2|) coincide. Notice that we have not computed the 
coupling Gp 1 p* n . In the infinite-mass limit it coincides with Gp**p n and therefore, in order 
to estimate the widths of the 1 + states, we assume that this equality holds for finite mass 
as well. From Eq. (|5.5| ) we obtain: 

-> D*n) ~ 165 MeV (5.6) 
T(B 1 -»■ 5V) ~ 360 Mel/ . (5.7) 

Also in this case we have taken mp 1 —mp* = 500 MeV (P = B, D) as suggested by HQET 
considerations. 

In order to perform some comparison with the experimental data, we write down also 
the formulae giving the widths of the = | + states flTTf : 

r(/-^^-) = s— x^W- («») 

6tt m P2 A x / w 

The strong coupling constant j£- can be estimated from the decay widths of the charmed 
state L> 2 (2460) -> Dtt,D*tt; using T tot (D 2 ) = 21 ± 5 MeV f28[, and assuming that only 



two body decays are relevant, one gets -jr- ~ 0.55 GeV 1 . From this result and from Eq. 
( p,10| ) one obtains for the state Df the total width T to t ~ 6 MeV; on the other hand 



the experimental width of the other narrow state observed in the charm sector, i.e. the 1 + 
£>i(2420) particle, is T tot (D 1 (2A20)) = 18±5 MeV |||. This discrepancy could be attributed 



to a mixing between the Di and the D[ states ||29|| . If a is the mixing angle, we have 



and therefore we get the estimate a ~ 16°. This determination agrees with the result of 



Kilian et al. in Ref. [11 



As for the B sector, evidence has been recently reported of a bunch of positive parity 
states B**, with an average mass mp** = 5732 ±5 ±20 MeV and an average width T(B**) = 



145 ± 28 MeV p0| , pl| . The observed states can be identified with the two doublets (2 + , 1 + ) 
and (l + ,0 + ). We note that the mass splitting, A = 500 MeV, between 5* and P states 
that we have chosen, agrees rather well with the experimental result in the B sector; our 
predictions for the widths, using the experimental value for the B** mass, are: T to t{Bo) — 
330 MeV, TtatiBj ~ 300 MeV, T tot (B 2 ) ~ 12 MeV and T^B^ ~ 10 MeV (we have 
neglected here the mixing which is a 1/rriQ effect). 

It is difficult to perform a detailed comparison of these results with the yet uncomplete 
experimental outcome; however, assuming that the result obtained by LEP collaborations 
represents an average of several states, its value is compatible with our estimate of the 
widths. Opal |3T] has also reported evidence of a B** state with mass m B ** = 5853 ± 15 
MeV and width r#»» = 47 ± 22 MeV. The width can be interpreted as connected to the 
decay B** — > BK, B*K. Assuming again that the width is saturated by two-particle final 
states, and using ms** = 5853 MeV, we obtain: 

T(B* s *(0 + )) ~ 280 MeV (5.12) 

1 + 

T(B* S * -> B*K) ~ 200 MeV (s t =- ) (5.13) 

3 + 

r(B**(l+)) ~ 0.45 MeV {s £ = - ) (5.14) 

r(S**(2+)) ~ 1.4 MeV . (5.15) 

Also in this case a detailed comparison with the experimental results cannot be performed 
without more precise measurements; we observe, however, that the computed widths of 
the different B** states are generally smaller than the corresponding quantities of the B** 
particles, a feature which is reproduced by the experiment. 

VI. EXCITED STATES CONTRIBUTION TO /d s //d 

Aim of this section is to study the contribution of the excited heavy mesons to the ratio 
of leptonic decay constants /d 8 //d- They are defined as 

< 0| dj^ 5 c\D+{p) >= if DVtx (6.1) 
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< 0\sj^ 5 c\D s (p) >= i/oJV (6.2) 

In the chiral SU(3) limit such a ratio is one, and the chiral corrections are expected to 
be of the order m s /A x . In terms of the fields H a and £ defined in Section II, the left-handed 
current q a l^ ~ ls)Q corresponds to [|3| 

L a ,= Z -F< 7m (1 - 75)^4 > + • • • (6.3) 

where the dots denote terms of higher order in the heavy mass and chiral expansion. At 
the lowest order, fjj a — fo — F / \jrrto- The effective couplings of the higher order terms, 
contributing to 5'L r (3)-violating corrections to the ratio /d s //d, are unknown. In addition, 
there are non-analytic corrections arising from chiral loops: in previous works p3| , |3"4]| the 



one-loop "log-enhanced" terms of the form m 2 log(m 2 //z 2 ) (m = m^, tuk or m v ) were kept, 
giving 

^ = 1 + 0.07 + 0.21/ (6.4) 

JD 

The effective coupling g, appearing in the heavy-light chiral lagrangian ( |2.5| ), gives the vertex 
D*Dtt, appearing in the loops. 

The loop corrections depend on an arbitrary renormalization point fi: this dependence 
is canceled by the //-dependence of the coefficients of higher order operators, which are here 
neglected. When fi is of the order of the chiral symmetry breaking scale A x «l GeV, these 
higher terms do not contain large logarithms and are supposed to be small compared with 
the ones coming from the chiral loops. 

The excited positive parity heavy mesons contribute to SU (3) violating effects as virtual 



intermediate states in chiral loops. In Ref. [12| the "log-enhanced" terms due to these 



excited-states loops has been computed: as we will see below, some of them are proportional 
to h 2 and others depend linearly on h. It has been pointed out that these terms could be 
numerically relevant and could invalidate the chiral estimate based only on the states D and 
D*; as we shall see below, however, the terms 0(h 2 ) and O(h), while important, tend to 
cancel. 
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In the following, we will present an independent calculation of the chiral loop contribu- 
tions to the ratio of leptonic decay constants, and we will give a numerical estimate obtained 
by using the QCD sum rules results for the couplings g and h. 

The chiral loop induced corrections to fo and fo s come from the diagrams of figs. 2,3 
and 4. 

The self-energy diagrams (fig. 2) give the following wave function renormalization factors: 



J D 



1 - 



3<£ 



3/2Ci(A D * D , A D * D ,m n ) + Ci(A D * D , A D * D , m K ) + -Ci(A D . D , A D * D , m v ) 



+ 



h 2 



[3/2C(A PoP , A PoP , mj + C(A Po3 p, A PosP , m K )+ 



+ -C(A Po p,Ap p,m v ) 



(6.5) 



Z Dr = 1 



3g 2 [ 2 
, r f, 9 2Ci(Ap;*p, s , A D * Ds ,m K ) + -Ci(A D * D , A D * D ,m v ) 



+ 



/) 2 



2C , (Ap )Ps , Ap oPs , m x ) + -C(Ap p, A Po p, 



(6.6) 



16tt 2 /2 

where the mass splittings Ap*p = Mp* — Mp, Ap*p s = Mp* — Mp s and Ap*p = Mp* — Mp 
are 0(l/rriQ), while the mass splittings Ap p = M Po — M Pl A PosP = M Poa — M P and 
Ap p s = Mp — Mp a between excited and ground states are finite in the limit ijiq — > oo. 
The functions C\ and C come from the loop integration and are defined as: 



(loOp = _>_ 

(2tt) 4 (q 2 - m 2 ){y ■ q - A)(v ■ q - A') ~ lQn 2 



(d(A, A', m)g aP + C 2 (A, A', m)v a v p ) 



and 



(6.7) 



C(A,A',m) = C l {A,A',m) + C 2 {A,A',m) 



(6.8) 



Performing the integration we obtain: 



d(A,A',m) 
C ( A,A',m) 



rrr 



9(A - A') 

2m 3 
9(A - A') 



A . A' , 

i?i( — ,m) - Hi{ — ,m) 
m m 

A A' 
H{ — ,m) — H{ — ,m) 

m m 



(6.9) 
(6.10) 
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where 



m 



Hi{x,m) = -12x + lOar 3 + (9x - Qx 6 ) log(— ) + 



12(x 2 - if 12 log(x + Vx 2 - 1) 
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,m 



H(x, m) = -9x A + (9x - -x) log(— ) + 



+ 18x 2 Vx 2 - llog(x + \lx 2 - 1) . 



(6.11) 



(6.12) 



The previous results are obtained in a renormalization scheme such that |+log(47r)— 7^+1 
0. 

For A = A', as for the wave-function renormalization factors, we find 



in 2 A 
C 1 (A,A,m) = — H[(-,m) 
9 m 

2m 2 A 
C(A, A, m) = ——H'(—, m) 
9 m 



(6.13) 
(6.14) 



where H'(x,m) = dH \^' m > . 



The diagram of fig. 3 is linear in h (the analogous proportional to g vanishes), and 
proportional to F + , defined in ( |3.29| ). 

Combining all the diagrams (fig. 2, 3 and 4 ) we obtain: 



D 



D 



32vr 2 / 2 



3 2 , / m ^\ 2 1 , m Kx 1 2, , m r)- 

-m w log(— ) + m K log(— ) + -m v log' 



/i 2 ' ' 6 " to > 2 



" //} ^Ci(A fl .j), A^d,??!,) + Ci(A D * D , A D * D ,m K ) + ^Cx(A D * D , A D * D , m v ) 



32tT 2 / 2 



+ 



+ 



h 2 D rs 



1 



32vr 2 / 2 
F + hp 

F 167T 2 / 2 

/d s = 



2 C(A Po p, Ap oP , m*) + C(Ap 0sP , Ap 0sP , m x ) + -C(A Po p, A PoP , m v ) 

3 1 
-C7(A Po p, 0, m^) + C(Ap 0sP , 0, m K ) + -C(A PoP , 0, m,,) 



(6.15) 



32vr 2 / 2 
32tt 2 / 2 

F 167T 2 / 2 



327T 2 / 2 



2<iog(^f) + ^io g (^|; 



.mz 



2C 1 (A D * Da ,A D * Da ,m K ) + -Ci(A£)»£), A^.^, m v ) 



2C(A PoPs ,Ap oPs ,m K ) + -C{A PoP ,A PoP ,m v ) 



2C(A PoPs , 0, m K ) + -C(A Po p, 0, m v ) 



(6.16) 
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From the previous formulae, using A PoP = 0.5 GeV, /i 
F = 0.30 GeV 3 / 2 , one gets numerically: 

F 



1, ^ 



Id 



F 



1 + 0.09 + 0.11/ - 0.33/i 2 - l.OO/i) 
, fl + 0.17 + 0.82/ - 0.66/i 2 - 1.15/1) . 



0.46 GeV" 3 / 2 and 

(6.17) 
(6.18) 



In the previous formulae we have kept only the leading order in the 1/uiq, i.e. we have put 
Ao»£) = in ( |6.15| ) and (|6.16|) : therefore we have to use for the couplings gr> and hp the 
asymptotic values g and h respectively. 

The value of g has been computed with QCD sum rules in [|17],|18|,[^] , giving a result in 
the range 0.2 — 0.4. From ( |6.17|) and ( |6.18|) , using g = 0.3 and h = —0.5, we get for the 
ratio of leptonic decay constants: 



L 

f. 



Ds 



1.09 . 



D 



(6.19) 



Without the contibution of the excited states, we would get /d s //d — 113: the contribution 
of the excited heavy mesons is slightly negative. Notice that the term in h 2 tends to cancel 
against the term linear in h; we also observe that its sign is unambiguously fixed by the sum 
rule (see Eq.( p.30|) , since the relevant quantity is the ratio hF + /F. 

The result ( |6.19|) is not very sensitive to the value of the mass splitting Ad d- For 
instance, if we take Ad d = 0.6 GeV we find joj fo = 1.12, while for A^^ = 0.4 GeV one 
obtains /d,//d = 1-05. We have also checked that (|6.19| ) depends weakly on the value of 
the renormalization scale \x. 

Keeping only the "log-enhanced" terms of the form m 2 log(m 2 //z 2 ), the ratio becomes: 



fDs 

fo 



32tt 2 /2 



,m 



k • 



.mz 



x M + 3g 2 + h 2 + 



1.06 



(6.20) 



The previous formula does not contain the parameter Ao d- The quoted number, 1.06, is 
for fi = 1 GeV: putting h = in (|6.20|), one obtains the results of P3J34J], i.e. Eq.QB.4p. A 
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different approximation has been put forward in [12[], where also terms A 2 log(A 2 //i 2 ) are 
kept, with numerical results similar to ours. 

In conclusion, the contribution of the excited heavy meson states to the ratio of leptonic 
decay constants fr>J fo is small and negative. The usual estimate of the £77(3) violation, 
including only the state D and D*, is not destabilized when including excited states, at least 
in this case. We stress that this result strongly depends on the sign of h: had it been positive, 
the final result would have been substantially different. Therefore we cannot exclude that 
for other observables the positive parity heavy mesons give a significant contribution to the 
chiral loops. In any case, we point out that the sign of the product F + /Fh, which enters in 
the formula (|6.20|) , is unambiguously determined by the sum rule. We also observe that the 



sign of h turns out to be negative also in ref. p5| , based on a chiral quark model. 



The numerical outcome Eq. fl6.20|) for fuj fn agrees with the theoretical results obtained 
by several groups by different models, e.g. QCD sum rules f36| , lattice QCD [[37]] and 



potential models fSBfl . As for the experimental results, we only have the upper bound |28] 

f D < 310 MeV (6.21) 
from MARK III collaboration [J39|, and the recent results from three experiments: 



f Da = 232 ± 45 ± 20 ± 48 MeV (WA75 [40]) 
f Ds = 344 ± 37 ± 52 ± 42 MeK (CLEO gTJ[) 

f Ds = A30t\lo ± 40 MeV ( BES El) • ( 6 - 22 ) 

Even though the theoretical value obtained for jo J jo is still consistent with data in 
( |6.21| ),( |BT2~2"|) , it should be noticed that the experimental results for fu s seem to indicate a 



value much larger than the theoretical estimates appeared in the literature p6|-p8|, which 
might signal a serious theoretical problem. In any event, better quality data are needed 
before any conclusion can be drawn. 
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VII. CONCLUSIONS 



We have computed the strong coupling constant of the positive parity heavy mesons 
Gb**bw and Gd**dtt, by QCD sum rules and light-cone sum rules. In the limit mq — > 
oo, the two methods give compatible results, but the 1/mg corrections are significant and 
unfortunately they are rather different in the two approaches. We have applied the results 
to the calculation of the hadronic widths of the positive parity B and D states: we have 
found that the calculated widths are in any case compatible with the recent preliminary 
LEP data on the orbitally excited B mesons. Furthermore, we have computed the chiral 
loop contributions of these states to the ratio /d s //d, and we have found that the chiral 
corrections consist of two sizeable quantities, which are however opposite in sign, so that 
the prediction for this ratio obtained using only the ground state heavy mesons is not 
significantly shifted. This cancellation is likely here to be fortuitous, and in view of the 
large value we have found for the coupling constant, leaves open the possibility that chiral 
contributions of the excited heavy mesons to other physical observables could instead be 
important. 
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FIGURES 




FIG. 1. The integration region of the double dispersive integral in Eq. ( p. 11 ). In the region 
D\ only the poles of the B, B* and B** particles are present, the region D2 includes the hadronic 
continuum. 
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FIG. 2. Self energy diagram 
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FIG. 3. Vertex correction involving positive parity heavy mesons 



FIG. 4. Vertex correction involving only light pseudoscalar mesons 
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